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Abstract. This paper deals with the higher dimension quasilinear parabolic-parabolic Keller-Segel 
system involving a source term of logistic type ut = V • (b(rt)Vu) — yV • (uVv) + g(u), rvt = A.v — v + u 
in Q x (0,T), subject to nonnegative initial data and homogeneous Neumann boundary condition, where 
0 is smooth and bounded domain in K", n > 2, <j> and g are smooth and positive functions satisfying 
ks p < (f> when s > So > 1, g(s) < as — gs 2 for s > 0 with <?(0) > 0 and constants a > 0, r, y, g, > 0. 
It was known that the model without the logistic source admits both bounded and unbounded solu¬ 
tions, identified via the critical exponent Y On the other hand, the model is just a critical case with 
the balance of logistic damping and aggregation effects, for which the property of solutions should be 
determined by the coefficients involved. In the present paper it is proved that there is 6q > 0 such 
that the problem admits global bounded classical solutions, regardless of the size of initial data and diffu¬ 
sion whenever ^ < 9$. This shows the substantial effect of the logistic source to the behavior of solutions. 
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1 Introduction 


In this paper, we consider the higher dimension quasilinear parabolic-parabolic Keller-Segel system with 
logistic source 


ut = V • (</>(u)Vw) — V • (ip(u)S7v) + g(u ), ( x , t) g 12 x (0, T), 


TVt = Av — U + It, 

<9?; 

<9n ’ 

u(x, 0) = u 0 (x), v(x,0) = v 0 (x), 


(x, t) € fl x (0, T), 

(x,t) € <912 x (0, T), 
ie!l, 


(i.i) 


where 12 C K" (n > 2) is a bounded domain with smooth boundary, r > 0. Functions (p, ip G C 2 ([0,oo)) 
satisfy 


if>(s) = xs (1-2) 

<p(s) >0, s > 0, ks p < </>(s), s > so, (1-3) 

with x > 0, k > 0, p e M, So > 1, and g G oo)) fullfills 

g(s) < as — g,s 2 , s > 0 (1.4) 


with <7(0)>0, and constants a > 0, g, > 0. Here, u and v represent the density of cells and the concentration 
of chemical signal respectively. The classical Kellcr-Segel system can be obtained by setting (f> = ip = 1 
and 5 = 0 in EH) , which models the mechanism of chemotaxis, and has been extensively studied since 
1970, we refer to 0[H (231C23 and the reference therein. 

Eq. ( 11 . 11 ) with g (it) = 0 is a type of refined models pursued by Hillen and Painter with the bacterial 
cells having a positive size, the so-called volume-filling effect. Beyond this, more general functions </> and V’ 
are involved to denote the diffusivity and chemotatic sensitivity, respectively mmm- When cp(s ) ~ s p 
and ip(s) ~ s q for large s, a critical exponent — on the interplay of (p and ip has been found to identify 
boundedness and unboundedness. Namely, if q — p < then all solutions are global and uniformly 
bounded [241 118) : however, if q — p > -4, unbounded solutions do exist [221; even finite-time blow-up may 
occur under some additional conditions n > 3 and q > i gcs]. 

Apart from the aforementioned system, a source of logistic type is included in (11.11) to describe 
the spontaneous growth of cells. The effect of preventing ultimate growth has been widely studied 
ng na Ea ng. In the related classical semilinear chemotaxis systems, that is when <p(u) = 1 and 
4 , { u ) — X u with x > 0, such proliferation mechanisms are known to prevent chemotactic collapse: In 
m ; for instance, it was proved that if r = 0 and g > - n 2 ^ + Xi solutions of the parabolic-elliptic system 
are global and remain bounded. The same conclusion is true for the fully parabolic system with r > 0 
if either n < 2, g > o m, or n > 3 and g > go with some constant go{x) > 0 [20j . This is in sharp 
contrast to the possibility of blow-up which is known to occur in such systems when 5 = 0 and n > 2 

EHHE21I22]- 

In this context, we intend to study ED with r > 0 under the conditions EH EH- It is our 
purpose to investigate the interaction among the triple of nonlinear diffusion, aggregation and the logistic 
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absorption. By taking u 1 1 as the test function to the first equation, and then substituting the second 
equation, the standard L 7 estimate argument yields 


1 d 
7 dt 


[ u 1 < -(7 - 1) [ 


M 7 +p - 2 |Vu | 2 + 


7-1 
7 + 9- 1 



— I —— [ u^- 1 Av 
7 + 9 - 1 Jn 



(1.5) 


Comparing the terms - / u 1+q and —p. / u 1+1 , it is easy to find that q = 1 is critical. It has 

7 + 9 — 1 Jn _ 

been proved that when q < 1 , the logistic dampening rules out the occurrence of blow-up regardless 
of diffusion [ 2 j. And when q > 1 , the strong diffusion with q— p < ensures global boundedness by 
[18] , without the help of the logistic damping. The critical case q = 1 is more involved: from m we 
may expect that under the balance of logistic damping and aggregation effects, the coefficients would 
determine weather the solution is bounded. In |3j, it has been proved that when q = 1 with r = 0, the 
solutions are bounded if /i > (1 — n (i_ p ) + )x f° r the parabolic-elliptic case. This makes an agreement 
with the above expectation. The main result of the present paper is the following theorem for the fully 
parabolic Keller-Segel system. 


Theorem 1. Suppose that O C 1" (n >2) is a bounded domain with smooth boundary, and Xit L i T > 0. 
Assume that ip(u) = x u > 4> an d 9 satisfy Then there is 9q > 0 such that if ^ < 9g, for any 

nonnegative ug £ C°(Cl) and vg £ W 1,r (fT) with r > n, Eq. El uniquely admits a classical solution 
{u,v) such that u £ C°(fl x [0,oo)) H C 2 ’ 1 (f7 x (0,oo)) and v £ C'°(f2 x [0,oo)) (~l C' 2 , 1 (f2 x (0,oo)) fl 
^ioc([0, ?max); W 1,r (Q)). Moreover, (u,v) is bounded infix (0, 00 ). 

Remark 1. We underline that the above result is independent of the value of p in (11. 3D . and thus extends 
the analogue result for the semilinear case ED]. Moreover, due to the technique used here, the convexity 
of f l (required in [ 20 ] ) is unnecessary in our theorem. 

Unlike using the trace embedding technique to estimate the boundary integral in mm , our ap¬ 
proach strongly relies on the Maximal Sobolev Regularity. 

The paper is arranged as follows. In section 2, we deal with the local existence and the extensibility 
of classical solution to El as well as a variation of Maximal Sobolev Regularity. Section 3 will be 
devoted to prove Theorem [1] 


2 Preliminaries 

The local solvability to El for sufficiently smooth initial data can be addressed by methods involving 
standard parabolic regularity theory in a suitable fixed point framework. In fact, one can thereby also 
derive a sufficient condition for extensibility of a given local-in-time solution. Details of the proof can be 
founded in [ 2 j. 

Lemma 2.1. Suppose f l C 1" (n > 3) is a bounded domain with smooth boundary, <f> and if satisfy 
E1"E3)> 9 fulfills J-M|), u 0 £ C°(fl) and vg £ W 1,r (fl) (with some r > n) both are nonnegative. Then 
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there exists (u,v) £ (C°(fi x [0, T max )) nC 2,1 (fl x (0,T max ))) 2 with T max £ (0, oo) classically solving (1.1) 
in ft x (0,T max ). Moreover, ifT mgx < oo, then 


limsup ||u(-,t)||i=(n) = oo. 

^/^max 


( 2 . 1 ) 


Note that in (EH), the term — - 


7-1 


. f m 7 +?-i Aw is unsigned. We thus intend to estimate its 
Jn 


_ 7 + <?- 1 

absolute value adequately. For this purpose, we will make use of the following property referred to as 
a variation of Maximal Sobolev Regularity, which will play an important role in the proof of our main 
result. The following Lemma is not the original version of a corresponding statement in [7, Theorem 3.1], 
but by means of a simple transformation by including an exponential weight function as in [lj. 

Lemma 2.2. Let r £ (1, oo), r > 0. Consider the following evolution equation 


rvt = Au — v + u, (x, t) £ fl x (0, T), 
dv 


dv 


= 0, 


(x, t) £ dVt x (0, T), 


( 2 . 2 ) 


j(x, 0) = xq(x), x £ Q. 


For each vq £ W 2 ’ r (Ll) (r > n) with = 0 on dLl and any u £ Z/((0, T); L r (Ll)), there exists a unique 
solution 

v £ W 1 - r ((0,T);L r (n)) nL r ((0,T);W 2 ’ r (ft)). 

Moreover, there exists C r > 0, such that if sq £ [0,T) ; v(-,Sq) £ W 2,r (fl)(r > n ) with = 0, then 

[ T [ e^|A v\ r <C r f [ e^ + are^o(|| v( . jSo) || Lr(n) + ||Ax(^o)ir L 7 n))- (2-3) 

J so J J so J ' ' 

Proof. Let w(x,s) = e s v(x,rs). We derive that w satisfies 


w s (x, s ) = Aw(x, s) + e s u(x, ts), (x, s) £ ft x (0, T), 
^ = 0, (x, s) £ dCl x (0,T), 

w(x, 0) = xq(x), x £ Cl. 


(2.4) 


Applying the Maximal Sobolev Regularity m Theorem 3.1]) to w, we obtain that 


[ [ \Aw(x,s)\ r <C r [ [ |e s u(x,rs)r+ C' r ||x 0 ||lr (n) +C' r ||Au 0 || 
Jo Jn Jo Jn 


L r ( n)- 


(2.5) 


Substituting v into the above inequality and changing the variables imply 


[ I e' s \Av\ r <C r f I e^ s w r + C' I .r||-y 0 |lL 7 n) + ^xllA-coll 
Jo Jn Jo Jn 

Consequently, for any sq > 0, replacing v(t) by v(t + sq), we prove (12.31) . 


zx( n)- 


□ 
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3 Proof of Theorem [T| 


In this section, we are going to prove our main result. Since the regularity obtained in (12.51) requires that 
the initial data satisfy homogeneous Neumann boundary conditions, we will perform a small time shift 
and thus use any positive time as the “initial time” to guarantee that the respective boundary condition 
is satisfied naturally. 

Specifically, given any s o G (0,T max ) such that so < 1, from the regularity principle asserted by 
Lemma [7TT1 we know that (u(-, So), v(-, s o)) G C 2 (fl) with = 0 on cftl, so that in particular we can 

pick M > 0 such that 


sup |K-, 0 )|| L =o (f2) < M, sup |M-,0)||i~ ( n) < M, and || Ac(-, so)|U~(fi) < M. (3.1) 

O<0<s o O<0<so 

Now we proceed to derive an a priori estimate which will constitute the main part of the work. 

Lemma 3.1. Suppose C K", n > 3, is a bounded domain with smooth boundary, r > 0 and x G M. 
For any 7 > 1, rj > 0, there exist p ltV > 0 and C = C( 7 , |f 2 |, p, y, 77, uq, vq) > 0 such that if p > p-y^, 
then 

IH-^llz/qn) < C 

for all t G (s 0 ,T max ). 


Proof. We fix so G (0, T) such that So < 1. For arbitrary 7 > 1, take u 7 1 as a test function for the first 
equation in ( 11 . 11 ) and integrate by part to obtain 


Id 


7 dt , 


(7 — 1) f u 1 2 (f(u)\Vu\ 2 + x(7 — 1) [ u 1 1 Vu • Vv + a f u 1 — p f u y+1 
Jo, Jo, Jo, Jo, 

- -- I Vu 7 • Vv + a f u 1 — p f u 1+1 

7 Jn Jn Jn 

X~ - [ « 7 A« + a f u 1 — p f u 7+1 

7 Jn Jn Jn 


7+1 
^7 Jn 


t 7 1 f 7 a if ,7+1 

u r — x - / u’/\v+[a + 


7 Jn 


T 7 J Jn 


u' — p / u 


7+1 


for all t G (so, Lmax). Here by Young’s inequality, for any e > 0, there exists ci > 0 such that 

(a + ^ + [ u 1 < e [ u 7+1 + ci(a,e, 7 )|H|, 

V +7 J Jn Jn 

where ci(a, £, 7 ) = A(1 + i) - ( 7 +i) £ - 7 ( a + ; AA) 7+1 . Young’s inequality also implies that 

— x ~—- [ « 7 A« <x[ it 7 |Ac| 

7 Jn Jn 


<77 / u 1+l + c 2 7? 7 X 


-7t,7+1 


|Au| 


7+1 


with C2 =supi(l +i)-( 7+1 ) < 00 . By substituting (13.31) and (13.41) into (13. 21) . we find that 


7>1 


d (\ 


dt \ 7 


7 <- 


7 + 1/1 


^ J u 7 ^ — (p — e — Tj) J u 1+1 + c 2 7 ? 7 X 7+1 


I Ac | 


7+1 


(3.2) 


(3.3) 


(3.4) 
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■ci(a,e, 7)|fi| 


(3.5) 


holds for all t G (sojT max ). Applying the variation-of-constants formula to the above inequality shows 
that 

1 r 1 r rt 
7 


■ f w 7 (-,t)<e ^ s °) — f u 1 (-,so) — (n — e — rj) I e 7 + ^ s ) f m 1+1 

Jq T Jq, jSn Jq 


+ C2V 7 X 7+1 


< — (p — e — rj)e 


-2±i(t-a) 


|Au | 7+1 + ci |£~2[ 


3 -^±i (t-,) 




SO 




+ C2V X e 


[ [ e^ ls |Au | 7+1 +c 3 (a,e, 7 , |fi|,s 0 ) 
J Sq ” 


for all t G (so,T max ), where 


C 3 (a, e, 7 , |fi|) = ci |Q| [ e 7 ^ 1{t s) + - / m 7 (-,s 0 ) 

J Sn ^ JQ 


is independent of t. Next, we apply Lemma [2721 to see that there is C 7 > 0 such that 


C 2 V 7 X 7+1 e * 


[ [ e^ s |Au | 7 + 1 

J S o « 

< c 2 ?7 -7 X 7+1 e _X ^‘ j^e^u^ 1 + ^re 2 ^ 80 ||u(-,s 0 )||^ 7+1(n) ^) 

Inserting m into m with some rearrangement, we finally arrive at 


1 

7 J o 


-s„,7+l 


C 3 


s J n 


■ [ M 7 (-, t)<-(fi-E-T)- c 2 C 1+1 pr 1 7 x 7+1 ) e A 
J n 

+ C 2 C 1 + 1 T V -^ X J+ 1 e- 2 ^ (t -^ |Kso)||&U (n) 

for all f G (so, Tjnax). Let /x 7j?? = rj + C 2 C 7 +if 7 _ 7 x 7+1 , we can choose EG (0,/z — /z 7iJJ ) such that 

p — £ — p — C2C 7 ?7 - 7 x 7+1 > 0 . 


(3.6) 


(3.7) 


(3.8) 


It is entailed that 

— [ u 7 (-, t) < C 4 (3-9) 

7 

for all t G (sq,T max ) with C 4 = c 3 + C 2 C' 7 T 7 _ 7 x 7 + 1 ||t'(so)liy 2 1 - T +i(Q)- This completes the proof by the 
above inequality together with m ■ □ 

Next, we invoke the well established Moser iteration to get boundedness of (u,v). 

Proof of TheoremUl By Morse’s iteration (Lemma A.l in [18]). we claim that there is 70 (n,p) > n > 0, 
determined via (A.8)-(A.10) in Lemma A.l of [IS], such that if 

IK-,*)IU 7 (n) < 00 (3.10) 
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for all 7 > 70 and all t £ (s o, T max ), then there exists C\ > 0 such that 

IK-,i)IU°°(n) < Ci (3-11) 

for all t £ (so, T max ). Actually, (Id.Kill implies that Vu is bounded. It is easy to check that all assumptions 
of Lemma A.l arc fulfilled. 

Let 9 0 satisfy 

hif u- 7,70 = inf (j? + c 2 C' 70+ ir/ _70 x 70+1 ) = ^~X- 

r}> 0 r/> 0 [70 

We see that ^ < 6q implies /r > with some 77 > 0. We know By Lemma GQ] that (Id. 101) holds, and 

hence (Id. Ill) is true. Combining with (Id. II) . we get that u is bounded in (0, T max ). The boundedness of v 
can be obtained by the standard parabolic regularity. Finally, Lemma [2.11 yields that (u,v) is global by 
contradiction. □ 
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